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Introduction 


One of the primary motivations behind this work is the attempt to understand the 
physics of a two-phase jet pump which constitutes part of a flow boiling test facility at 
NASA-Marshall. The flow boiling apparatus is intended to provide data necessary to 
design highly efficient two-phase thermal control systems for aerospace applications. 
The facility will also be capable of testing alternative refrigerants and evaluate their 
performance using various heat exchangers with enhanced surfaces. The test facility is 
also intended for use in evaluating single-phase performance of systems currently using 
CFC refrigerants. 

Literature dealing with jet pumps is abundant and covers a very wide array of 
application areas. Example application areas includes vacuum pumps which are used in 
the food industry, power station work, and the chemical industry; ejector systems which 
have applications in the aircraft industry as cabin ventilators and for purposes of jet thrust 
augmentation; jet pumps which are used in the oil industry for oil well pumping; and 
steam-jet ejector refrigeration, to just name a few. Examples of work relevant to this 
investigation includes those of Fairuzov and Bredikhin (1995). 

While past researchers have been able to model the two-phase flow jet pump 
using the one-dimensional assumption with no shock waves and no phase change, there is 
no research known to the author apart from that of Anand (1992) who was able to account 
for condensation shocks. Thus, one of the objectives of this work is to model the 
dynamics of fluid interaction between a two-phase primary fluid and a subcooled liquid 
secondary fluid which is being injected employing atomizing spray injectors. The model 
developed accounts for phase transformations due to expansion, compression, and 
mixing. It aslo accounts for shock waves developing in the different parts of the jet pump 
as well as temperature and pressure dependencies of the fluid properties for both the 
primary two-phase mixture and the secondary subcooled liquid. 

Analysis 

The objective of this research effort is to develop an analytical model capable of 
predicting the performance of a two-phase jet pump. The model developed will account 
for phase transformations which are likely to take place in the nozzle, mixing chamber 
and diffuser (see Figure 1). Complete details of this model can be obtained from the 
NASA colleague. 

Primary Nozzle . The approach employed in the analysis described in this report is an 
extension of the isentropic homogeneous expansion (IHE) model described in the 
literature. In addition to the primary nozzle, the extension to the model covers the rest of 
the jet pump components and is thus labled as an isentropic homogeneous expansion/ 
compression (IHE/C) method. This method assumes that the velocities of the two phases 
are equal and that thermal equilibrium between the phases exists. It further assumes that 
both the expansion in the nozzle and compression in the diffuser are isentropic and that 
the working fluid property data correspond to those of a static, equilibrium, two-phase 
system with plane interfaces. 
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In order to be able to deal with the flow in the jet pump as a homogeneous two- 
phase mixture, equivalent properties are derived for the mixture in terms of 
corresponding properties for its vapor and liquid phases. In the nozzle part of the jet 
pump, the two-phase primary fluid will be treated as predominantly gaseous with liquid 
particles suspended in it. It is hypothesized that the presence of the liquid particles in the 
mixture is such that the latter can still behave as a gas, but with modified properties. 

The mass flow rate of the two-phase fluid flowing through the primary nozzle can 
be expressed according to the following equation (Shapiro, 1953): 
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where the quantities with the subscript “n” are those of the two-phase primary fluid while 
flowing through the nozzle. To obtain an expression for the specific heat ratio of the 
primary fluid, y n , the mixture quality in the primary nozzle, x n , is used in conjunction with 
the definition of y n : 
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To obtain an expression for the two-phase gas constant of the primary fluid, R n , we utilize 
the expressions for Cp n , C v , n and the definition of R„ to get 


With the exception of the specific heat at const int-pressure of the vapor phase, the 
properties of interest are fairly independent of temperature and pressure. For the two- 
phase primary fluid, the properties C p>n , C v ,„, Yn. *nd R n are therefore going to be 
predominantly functions of the primary fluid vapor-liquid composition or the quality, x n , 
and will be functions of temperature to the extent that C p , gjn depends on the temperature. 
Having established this fact, we can now proceed to perform two-phase calculations 
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throughout the nozzle to establish the conditions of the two-phase primary fluid at the 
nozzle exit. This will be done by first computing the mixture specific entropy at nozzle 
inlet, s n j, as follows: 

S ni = S ni.L + X n, ( S ni, g ~ \,.l) ( 7 ) 


where s n i, g and s n i,L represent the saturated vapor and saturated liquid specific entropies 
evaluated at the nozzle inlet pressure, p n j. 

Next, the nozzle exit pressure is assumed and the saturated vapor and saturated 
liquid specific entropies, s ne ,g and s ne ,L, respectively, are computed. Since the flow is 
assumed isentropic throughout the nozzle, the specific entropy at the nozzle exit will be 
equal to that at the nozzle inlet. This allows the quality of the primary fluid at the exit, 
Xne, to be computed 

( S ne~ S ne,L^ / - 0 \ 
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Since the exit pressure has been assumed, the exit temperature will have also been 
specified (since the fluid is in the mixture region). This allows computing the vapor 
specific heat at constant pressure, C p ,g, ne , the gas constant of the vapor phase, R g , ne , and 
the vapor phase ratio of specific heats, y g ne all evaluated at the nozzle exit. These, in 
addition to the primary fluid quality computed from Equation (8), can then be used to 
compute the primary fluid specific heat ratio from Equation (5). 

Now that the exit pressure and mixture specific heat ratio are known, the exit 
Mach number can be computed employing the following equation (Shapiro, 1953): 
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Using the exit Mach number, M ne , along with the value of the primary fluid specific heat 
ratio, y ne , the area ratio between the exit and throat can be computed using the following 
equation (Shapiro, 1953): 
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Since A ne /A m is a known quantity for a specified nozzle geometry, the area ratio 
computed from Equation (10) can now be compared to that assigned by the nozzle 
configuration. If the two values did not agree, then Equation (10) can be used to compute 
the Mach number at the nozzle exit corresponding to the actual area ratio, but employing 
the computed value of Yne from the previous calculation step. The new value of the Mach 
number can now be used to compute a new value for the nozzle exit pressure, p ne , using 
the following equation (Shapiro, 1953): 
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Again, the value of y ne used in Equation (11) is the same as the one employed in the 
previous calculation step. Once the nozzle exit pressure has been found, the exit 
temperature, T ne can also be found since the primary fluid is a two-phase mixture. The 
new exit temperature is then used to correct the value of the vapor phase specific heat at 
constant pressure, Cp igine , and the vapor phase gas constant, R g , ne , both evaluated at the 
nozzle exit. As mentioned earlier, variations in the gas phase ratio of specific heats as 
well as the specific heat of the liquid phase are relatively insensitive to temperature 
variations in the temperature range of interest to this study. While updating the values of 
these parameters to account for the newly calculated temperature value at the nozzle exit 
is relatively straightforward, it will not be performed here since its effect on the accuracy 
of the results should be fairly minimal. 

With the newly calculated exit pressure, the specific entropies of the saturated 
vapor and saturated liquid evaluated at the nozzle exit pressure can be found and used to 
calculate the updated value of the primary fluid quality at the nozzle exit using Equation 
(8) and the isentropic flow assumption. This quality can, in turn, be used to calculate the 
primary fluid ratio of specific heats at the exit, y ne , employing Equation (5). Knowing, 
Y„e, a new value for the Mach number at the exit can be calculated using Equation (9) and 
subsequently a new area ratio can be calculated using Equation (10). The calculated area 
ratio should then be compared to the actual area ratio, and in case they do not agree, the 
whole procedure is repeated until convergence. 

After convergence has been reached, the values of all the relevant parameters of 
the primary fluid at the nozzle exit should be updated. This should include the ratio of 
specific heats y ne , gas constant R ne , specific heats C p , n e and C v . n e, temperature T ne , pressure 
p ne , quality x ne , entropy s ne , enthalpy h ne , density p ne , and Mach number M ne . The 
enthalpy can be computed using the following equation: 

Ke = KeX + X ne ( h ne, g ~ Ke.D ( 12 ) 


where h ne ,L and h ne ,g represent the saturated liquid and saturated vapor enthalpies 
evaluated at the exit pressure, p ne - The density, p ne , on the other hand, can be computed 
from the following equation (Shapiro, 1953): 
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Of course the specific volume at the nozzle exit, v ne , is the reciprocal of the density and 
can easily be determined. The mass flow rate of the primary fluid can now be computed 
from Equation (1). The area employed in the equation should correspond to the exit 
section since the Mach number and fluid properties used correspond to the exit. 
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With the knowledge of the mass flow rate, as well as the exit area and density, the 
nozzle exit velocity can be computed using the continuity equation 

V ne = mj (p ne A ne ) (14) 

where the quantities v n j >g and v n j_L represent the specific volume of the saturated vapor and 
saturated liquid of the primary fluid evaluated at the nozzle inlet pressure, p n i. 

Liquid Sprayers . As was reported earlier, atomizing sprayers are used to inject 
subcooled liquid into the mixing chamber of the jet pump. In order to relate the pressure 
and velocity of the secondary liquid at the inlet of the mixing chamber to those at the inlet 
of the sprayers, Bernoulli’s equation is used 



Additional details can be found with the NASA colleague. 

Mixing Chamber . The governing conservation equations for a control volume 
representing the mixing chamber are those of the continuity, momentum, and energy (in 
integrated forms): 

m m =m n + m s (17) 

K P ne + A, p mi cos a - A me p me = m m V me - m n V p mi - m s V s mi cos a (18) 




(19) 


Additional details about the mixing chamber analysis can be found with the NASA 
colleague. 

Throat and Diffuser . Because it is highly likely that a condensation shock would take 
place some where in either the throat or the diffuser, the latter will be divided into three 
segments; the part of the diffuser upstream of the shock, the shock layer itself, and the 
part of the diffuser downstream of the shock. For both parts upstream and downstream of 
the shock, control volume analyses will be performed. For the shock layer, the standard 
normal shock relations for a two-phase substance will be implemented. 

The problem on hand can proceed in two different but equivalent routes. The first 
route is to assume the location of the condensation shock and let the solution predict the 
diffuser exit conditions. This route leads to predicting the jet pump compression ratio 
with a prescribed shock location. The other route is one in which the exit pressure of the 
diffuser is controlled (such as the work of Fairuzov and Bredikhin, 1995) and the location 
of the shock is determined as part of the solution. In this report, the latter approach is 
chosen. Complete details on the governing equations of the throat and diffuser can be 
found with the NASA colleague. 


XLI-5 



Conclusions 


The research effort on which this document partly reports described a relatively 
simple model capable of describing the performance of a two-phase flow jet pump. The 
model is based on the isentropic homogeneous expansion/compression hypothesis and is 
capable of fully incorporating the effects of shocks in both the mixing chamber and the 
throat/diffuser parts of the pump. The physical system chosen is identical to that 
experimentally tested by Fairuzov and Bredikhin (1995) and should therefore relatively 
easy to validate. 
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ni » nozzle inlet 
nt - nozzle throat 
ne ■ nozzle exit 


mi ■ mixing chamber inlet x ■ upstream of shock 

me ■ mixing chamber exit y - downstream of shock 

di - diffuser inlet 

de - diffuser exit 


Figure 1. Physical model of the jet pump 
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